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Abstract 

We establish the relation between two objects: an integrable system related to 
Painleve II equation, and the symplectic invariants of a certain plane curve Htw- This 
curve describes the average eigenvalue density of a random hermitian matrix spectrum 
near a hard edge (a bound for its maximal eigenvalue). This shows that the s — > — oo 
asymptotic expansion of Tracy- Widow law Fque(s), governing the distribution of the 
maximal eigenvalue in hermitian random matrices, is given by symplectic invariants. 
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Introduction 



The Tracy- Widom law Fque(s) governs the distribution of the maximal eigenvalue in 
large hermitian matrices, sampled randomly in the GUE ensemble [32]. In fact, this 
law is universal, it is also valid (under some hypothesis) for non gaussian ensembles 
in the unitary symmetry class [9], for large Wigner matrices [30] and for a certain 
class of large random complex covariance matrices [31]. Just like the Gumbel, Frechet 
and Weibull distribution are the possible universality class for the maximum of a large 
number of independent variables [16, 20], Fgue is a possible universality class for the 
maximum of strongly correlated variables. Thus, it is a very important law in physics 
and mathematics, and it has been observed in numerous phenomena and experiments 
[24, 3, 29]. Tracy and Widom have written in 2002 a concise review [33] and provide 
article references for the situations where Fque(s) appears. 

The GUE ensemble [25] is defined as the space of N x N hermitian matrices Hn 
endowed with probability law: 

d(i N (M) = idMe-"^ 
Zn 

where dM is the canonical Lebesgue measure on the real vector space Hn- The decom- 
position M = U diag(Ai, . . . , \^) W induces the probability law on the eigenvalues: 

1 N A 2 

d/i7v(A) = — — - \Xi - Xj\ 2 • TTdAi e~ N ^ (0-1) 

l<i<f<N ti 

The expected maximum eigenvalue is limjv-^oo E(A max ) = 2y/i, and the scale of fluctu- 
ations of the maximum around this value is N~ 2 ^ 3 . Here, we set t — 1. The following 
limit exists and defines the GUE Tracy- Widom law F G ue : 

F GUE (s) = lim Prob[A max < 2 + s N" 2 ^] (0-2) 

In their original paper [32], Tracy and Widom were the first to characterize this law: 

Theorem 0.1 Let q(s) be the Hastings-McLeod solution of Painleve II, uniquely de- 
fined by: 

q" = 2q 3 + sq, q{s) ~ \f^j2- 

s— >— oo 

If we define H(s) = f s °° q 2 (o~) da, 



Fgue (s) = exp (- £ H(a) da^ . 
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H(s) can be identified with a Hamiltonian for PII [27], and Fgue(s) with a r-function 
associated to this family of Hamiltonians ([18] in the sense of Okamoto, [5] in the sense 
of Jimbo-Miwa-Ueno). From Theorem 0.1, Fque(s) can be shown to have the following 
asymptotic expansion when s->-oo: 

Fgue(s) = C exp ( - kj£l + ]T(- s /2) 3 ^ AA . (0-3) 

In some sense, A = — |, A x = — |, and the higher order coefficients A g can be 
computed recursively thanks to Painleve II equation. This method does not give the 
constant C = 2 1 / 24 e^^~ x \ which was obtained later by Deift et al. by Riemann-Hilbert 
asymptotic analysis [10]. This article aims to provide an alternative description of A g : 

Proposition 0.1 A g are the symplectic invariants F ff (E PII ) associated to the spectral 
curve of equation: 

(E PII ) : y 2 = ±x 2 (x 2 + 4). 

We will explain briefly the notion of spectral curve S and symplectic invariants F 3 (S) 
(Section 1.1). Those algebro-geometric objects were introduced axiomatically in [12]. 
One of their important property is their invariance under transformations of the spec- 
tral curve (x, y) — > (x , y Q ) such that \dx Q A dy \ = \dx A dy\. So, we have equivalently: 



Proposition 0.2 A g are the symplectic invariants F 9 (S TW ) associated to the spectral 
curve of equation: 

(S TW ) : y 2 = x+--2 

x 

In the chronology, we first obtained Proposition 0.2 heuristically by studying large 
deviations in an ensemble of random hermitian matrices conditioned by A max < a [6]. 
Then, we realized that the curve Etvk is equivalent up to symplectic transformation 
to Spu, and the latter appears in relation with Painleve II [17]. Now, we give a direct 
proof of Proposition 0.1 relying on earlier results in integrable systems. 

Meanwhile, we felt the need to prove the claim that Tau functions for integrable 
systems are related to symplectic invariants. As far as Tracy- Widom law is concerned, 
we only need this claim for 2x2 Lax pairs, and it was proved in a longer version of 
this article [7]. We are currently working on the generalization to d x d Lax pairs. 
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Outline 



We begin with an introduction to "symplectic invariants" and the formalism of the 
"topological recursion" (Section 1.1). We explain its relation to integrable systems 
in the case of 2 x 2 systems (Section 1.2-1.3, with Thm. 1.1 and 1.2 as main tools). 
Next, we show how this can be applied to an 2 x 2 integrable system whose compati- 
bility condition is the Painleve II equation, and which admits Fgue(s) as Tau function 
(Section 2). This will prove Proposition 0.1. 
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1 Topological recursion and 2x2 integrable systems 



1.1 Axiomatic of the topological recursion 

Let us recall briefly what the topological recursion consists in (for a extensive review, 
see [14]). It is an algorithm associating some numbers F 9 (E) and differential forms 
oj9(E)(zi, . . . , z n ) to a regular spectral curve E. For our purposes, we call spectral 
curve the data of: 

• a plane curve (C,x,y), in other words a Riemann surface C, with y and dx 
meromorphic functions C — > CP 1 (In most of our examples, C = CP 1 is the 
Riemann sphere, and thus x and y are complex rational functions). 

• a maximal open domain U C C on which x is a coordinate patch, called physical 
sheet. 

• a Bergman kernel P>(zi, z 2 ), i.e. a differential form in z\ G C and in z 2 6 C, such 
that, in a local coordinate £: 

= ,5r )d ft +ow 

*i-«a ({(zx) - £(z 2 )) 2 

The zeroes of dx lying on dU are called branchpoints (name them a^), and a spectral 
curve is said to be regular when these zeroes are simple and dy(a,i) 7^ 0. In other words, 
when yfx — x{ai) is a good coordinate around a ; 6 C and y — y(aj) oc Wx — x(cti). 
Before giving the full definitions of w 9 (E) and F 9 , we need to introduce: 

• The local conjugation. For 2 in a neighborhood of a;, there exists a unique 2^ z 
such that x(z) = x(z). 

• The recursion kernel. We define, for zq G C, and z in a neighborhood of a;: 

(2/(2) -2/U))dx(z) 
i? is a differential form in z 0; an d the inverse of a differential form in z. 

• For 2 in a neighborhood of dj, we introduce an (arbitrary) primitive of ydx. 
Then, we construct: 

= -y(z)dx(z), w^zx, * 2 )(E) = z 2 ) (1-2) 
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and by recursion on M = 2g + n — 2 > 0: 



j9(zi,z 2 ,...,z„) = V? Res _R(zi,£) 
«■ ' ^-^f^ai 



KC,l,0<Kj 

(1-3 



means that we exclude the terms where cu^ itself appears. Eventually, for g > 2: 

A — Zq z — ' 

We refer to [12] for the construction of F 1 and F°, which is more involved. Let us state 
the main properties of F 9 (S) and cj^(E): 

• For 2 - 2# - n < 0, w»(E) G T*(C) ® ■ ■ ■ ® T*(C), i.e. . . . , z n ) is a 
meromorphic differential form in each Zi G C, symmetric in all z^s, and it has 
poles only at the branchpoints. There is no residue at these poles, and their 
maximal order is 2(3g + n — 2). 

• F 9 (E) is invariant under any transformation (x, y) — > (x Q , y ) such that |dxAdy| = 
|dx D A dy a \. For this reason, the F 5 (S) are called symplectic invariants 3 . 

• cu^(S) have nice scaling properties when the spectral curve approaches a singular 
spectral curve. We will be more precise when needed. 

• c<^(£) have nice properties under variation of the spectral curve. 
1.2 Relation to 2 x 2 integrable systems 

Let L(x) be a 2 x 2 matrix whose entries are rational in x. We consider a solution ty(x) 
to the differential system: 

±d x *(x) = L(xMx), *=(t I) (1-5) 

N is a parameter of L. The 1/N in front of the derivative is here for convenience, and 
can always be absorbed in a redefinition of x and L. It is not restrictive to assume 
TrL = 0. So, d x (detV) = 0, and we can choose the normalization: 

det * = 1 



3 However, cj^(S)(zi, . . . , z n ) is shifted by an exact differential form in each Zj under symplectic 
transformations 
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To this system is associated the Christoffel-Darboux kernel: 

1p(Xi)4>(x 2 ) - ljj(x 1 )(f)(x2) , , 

lC(xi,x 2 ) = (1-6) 

X\ - x 2 

In [4] were also introduced correlators W n (ii, . . . , x n ) and connected correlators 
W n (xi, . . . , x n ). The connected correlators are defined by: 

Wi(x) = lim f ]C{x,x') — 



W 2 (xi,x 2 ) = -JC(xi,x 2 ))C(x 2 ,xi] 



(xi - x 2 ) 2 



a cycles of & n i=l 

and the correlators by: 

W n (xi,... J x n ) = "det"/C(x ij z J -) (1-8) 

where "det" means that each occurence of JC(xi,Xi) in the determinant should be 
replaced by Wi(xj), and each occurence of K,{x ir Xj))C(xj, x^) by — W 2 (xi, Xj). For 
example: 

Wi = V 9*0 -il)d x (f) = - (d x ip 4> - d x ip 4>) 

It can be checked that the correlators are symmetric in the Xi's, and that they do not 
have poles at coinciding points x% — Xj, i ^ j. The spectral curve of a first order 
differential system is defined by the plane curve of equation: 

det(y - L(x)) = (1-9) 

(since this is an algebraic equation, this defines a Riemann surface C and two mero- 
morphic functions on it, related by the algebraic equation). As physical sheet, we take 
the open maximal domain U C C containing one preimage of the poles of L(x) (here 
x = oo), and as Bergman kernel: 

B{zx, z 2 ) = —dx(z 1 )dx(z 2 ) K,{x{zi), x(z 2 )) K,(x(z 2 ),x(z 1 )) (1-10) 

These objects can be studied with the following theorem. 

Theorem 1.1 Assume that: 

(i) L has a limit when N — > oo. 
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(ii) The spectral curve Y>n of the system Eqn. 1-5 has a large N limit Eoo which is 
regular, and has genus 4 0. 

(Hi)' Wi(x) admits an asymptotic expansion when N — > oo of the form : W± = 
^ g>0 iV 1_2s Wf , and for g > 1, each coefficient Wf(x) for g ^ may have 
singularities only at branchpoints of E^ . 

Then, W n (xi, . . . , x n ) admits an expansion of the form: 

w n = J2 N2 ~ 29 ~ n w° (1-11) 

The expansion coefficients of the correlators have only singularities at the branchpoints 
o/Eqo, and are computed by the topological recursion applied to E^: 

W^(x(z 1 ),...,x(z n ))dx(z 1 )---dx(z n ) = ^(E^)^, . . . , z n ) 

dxpi) dx(z 2 ) , . 
~ 5 ^ 2 (x( Zl )-x(z 2 W ( } 

The result was first proved in [4] with a stronger hypothesis (Hi), which was weakened 
to hypothesis (Hi)' in [7]. Let us explain its origin in brief. Following [1], we can always 
embed d x ^/ = L\l/ in a hierarchy of compatible differential systems, with a full family 
of times t = (tj)j e ^: 

^d x V(x,?)=L(x,i)y(x,i), ^d t .y(x,i) = M j (x,i)V(x,i) (1-13) 

Then, one can show that the expansion of Eqn. 1-11 and the analytical structure of its 
coefficients for all n > 1 is implied by (Hi). The correlators W n were defined such that 
the coefficients of the expansion satisfy loop equations. Actually, the solution of 
loop equations with prescribed analytical structure is unique. This analytical structure 
is encoded in the large N limit spectral curve E^. And, the topological recursion was 
designed precisely to produce the unique solution to the loop equations with analytical 
structure given by Eqo. 

1.3 Tau function and symplectic invariants 

It was claimed in [4] and justified in [7] that the Tau function of the integrable system 
Eqn. 1-13 is given by a resummation of the symplectic invariants of Eqo. 

Theorem 1.2 We make the following assumptions: 



4 A similar result holds when Soo is not of genus under an extra hypothesis on Wf . This will not 
be needed here. 
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[i-iii') The hypothesis of Thru. 1.1. 

(iv) Eqn. 1-13 has a solution ^/(x,i) which behaves like: 

t) ~ (1 + o(l))e NTaix d (1-14) 
in a given sector when pole of h(x,F). 

(v) T a (x,t) is a 2 x 2 diagonal matrix. 

Then, the Jimbo-Miwa-Ueno Tau function, r(t), associated to that solution admits has 
the same large N asymptotic as: 

Fffj = C exp (j2 N2 ~ 29F9 P°°&n t 1 " 15 ) 

\g>0 / 

up to a constant C which does not depend on t. 

2 Application to Painleve II 

2.1 Integrable system associated to Painleve II 

The Painleve II equation q"(s) = 2q 3 (s) + sq(s) appears [19] as the compatibility 
condition of the following system for s): 

= 

= (2-1) 
= \l> exp \i (|x 3 + sx) o~3~\ when x — > +oo 

where $f — 1 + 0(l/x) when x — > +oo. The Lax pair (L, M) is given by: 

, , f -4tx 2 - i(s + 2q\s)) Axq(s) + 2ip{s) \ , , 

1 ' ' V 4xq(s) - 2ip(s) Aix 2 + i{s + 2q\s)) ) 1 > 

M(x,s) = q ( s) ) (2-3) 




q(s) ix 

The necessary condition of existence of \1/ is <9 S L — d x Nl + [L, M] = 0. This implies 
that q(s) is solution of Painleve II, and p(s) = q'(s). The asymptotic behavior of \I/ 
(Eqn. 2-1) determined the asymptotic behavior of q{s), namely: 

Q{s) ~ Ai(,) ~ ^pf/P (2-4) 

Hastings and McLeod [21] have shown that requiring this asymptotic determines a 
unique solution of Painleve II. Moreover, this solution is also the unique one with left 
tail asymptotic: 

q(s) ~ y^I/2 (2-5) 

s— oo 
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The existence of the solution ^f(x,s) was shown in [19]. This system is an example 
of isomonodromy problem, itself closely related to a Riemann-Hilbert problem. Many 
authors have contributed to the study of these systems, and we refer to [17] for a review 
of the theory. In this section, our goal is to do some simple checks to put ourselves in 
the framework on Section 1. 

Let us introduce a redundant parameter N. We define: 

{ iZtffls ' 9{s) = N^QiS) (2-6) 
Then, Eqn. 2-1 is equivalent to: 
±d x V = 

±d s V = (2-7) 
* = * exp [IN (§A 3 + SX) a 3 ] when X -> +oo 



with the Lax pair 
L(X,S) 



-AiX 2 -i{S + 2Q 2 {S)) AXQ(S) + 2 ^^ 

2iQ'{S) A - v2 ,-f Q ,0/02/ 



AXQ(S)-^^ AiX 2 + i(S + 2Q 2 (S)) 

M < A ' S > - (qIs) V) < 2 - 8 > 

and the compatibility equation: 

j^Q"(S) = 2Q(S) 3 + SQ(S) (2-9) 

Again, Q(S) is given by the Hastings-McLeod solution of Painleve II, which is such 
that (see Eqn. 2-5) \imN^ 00 Q(S) 2 = — |- for any S < 0. The reason for introducing N 
is that the s — > — oo asymptotics for the system Eqn. 2-1 corresponds to the N — > oo 
asymptotics for the system of Eqn. 2-7. We have now the notations of Thm. 1.1. 

2.2 Spectral curve 

After Eqn. 1-9, the finite N spectral curve for this system is: 

Y 2 = -16A 4 - 8X 2 S -S 2 - AQ\S) - ASQ 2 (S) + 4 ^if^ (2-10) 

N z 

In the large N limit, since Q 2 (S) ~ —S/2 (we assume S < 0), we obtain: 

(Boo): Y 2 = -16A 2 (x 2 + £) (2-11) 
It can be brought in a canonical form by rescaling: 

(£ P „) : y 2 = ^x 2 {x 2 + 4) (2-12) 
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where 

x=z \/5 x ' y=l ~s Y (2_i3) 

For this transformation, YdX = (—S/2) 3j/2 ydx. Going from the topological recursion 
of Eqo to that of E PII is only a matter of rescaling: 

F^Eoo) = (S/2) 3 ^ F*(E P „) (2-14) 

4 ff) (^oo) = (-5/2) 3 ( 1 "^/ 2 )^(E PII ) (2-15) 

Thus: 

= ^(iy 2 / 3 (-S/2)) 3(1 " 9) ^(E PII ) 

= £(- S /2) 3M F*(£ P „) (2-16) 
2.3 Proof of Proposition 0.1 

In our case, L(X, S) has only one pole, at X = oo, which is already present in 
L/°)(X, S). We define Y(x) by Eqn. 2-11 with the choice of the branch of the square 
root imposed by the X — > oo asymptotics in Eqn. 2-7: 



Y(X) = -MX \/X 2 + -, Y(x) ~ -AiX 2 (2-17) 
V 2 x— >+oo 

If the hypothesis of Theorem 1.2 were satisfied, it would prove the identity (in the sense 
of s — > — oo asymptotic) between 

lnF GUE (s) = lnr(s) and ^(E^) = ^(- s /2) 3 ^) F*(E PII ) 

which is the content of Proposition 0.1. For the moment, the state of our checklist is: 

• Hypothesis (z), (iv) and (v) are automatic with Eqn. 2-7. 

• Eqo is a regular spectral curve of genus 0, i.e fulfils hypothesis (ii). This is readily 
seen with the rational parametrization: 



Wi has an expansion in odd powers of N (see Lemma 2.2 below), and Wf has 
only singularities at branchpoints of Eqo (see Lemma 2.1 below). Hence (Hi)'. 
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We now turn to the technical lemmas. 



Lemma 2.1 if)(X, S) (resp. ip(X,S), 4>(X,S), 4>(X,S)) admits a 1/N expansion: 

'4 



^(X, S) = ^ \X, S) + J2 S) exp 



i>i 



( ^X 3 + sx 



and for all I > 1, tp^(X, S) has only poles at X = ±y (the branchpoints of Tjqo). 
In particular, it does not have singularities at X = 0, the other zero ofY. 

proof: 

We only sketch the proof. We already know that Q^°\S) = \J~^§-, and the Hastings- 
McLeod solution of Eqn. 2-9 has a 1/N 2 expansion: 

Q(s) = q {0 \s) + J2n- 21 q {l \s) 

1>1 

Hence, L and M have a 1/N expansion. The large N limit spectral curve associated 
to the system dx^ = is: 



Y(X, S) = ±AiX^X 2 + S/2 
and -0 to leading order is given by: 



tp (0) (X,S) = cte (0) (5) 
= cte^(S) 



2Y(X) 
-S/2 



exp 



X 2 + S/2 



-N I F(£)d£ 
1 



exp 



TiN ( ^X 3 + SX 



To agree with the asymptotic required in Eqn. 2-7, we must choose the minus sign. 
On the other hand, the large N limit spectral curve associated to the system 

i d s V = Mtf is: 

Y(X, S) = iV-detMW = ±i^X 2 + S/2 

Hence, ip is also given to leading order by: 



N / Y(X,a)da 



^(X,g)=cteW(X) ^|^ ) exp 

One can determine recursively the solution of each differential system (with respect to 
X, or with respect to S), and one finds an expansion of the form of Eqn. 2-18. By 
plugging this expansion alternatively in the two systems, it is possible to show that 
can have singularities only at common zeroes of Y(X, S) and Y_(X, S). In other 
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words, the possible singularity at X = obtained with the system with respect to X, 
is shown to disappear if we use the system with respect to S. At all orders in 1/N, 
only the singularities at the simple branchpoints of remains, i.e. at X = ±^S/2. 
The discussion for 0, and ip is similar. □ 

Lemma 2.2 Wi(X) admits an expansion in odd powers of 1/N: 

Wx(X) = J2 Nl ~ 2g Wf( x ) 

and the singularities ofWf(X) for g ^ are found only at branchpoints X = ^J—S/2. 
proof: 

From the asymptotics of ip, we see that Wi has at least an expansion in powers of 1/N, 
and the position of the singularities away from X = oo at all orders is a consequence 
of Lemma 2.1. One can also check that a singularity at X = oo can only appear in 
the leading order W°. Let us stress the dependence in iV by writing L^r, and for 
the solution of Eqn. 2-7. We observe that *L_tv = L^v, which implies that *\l/~^ v obey 
the same differential system as Moreover, '^1^ has the same asymptotic behavior 
near x — > oo as ^n, and is also of determinant 1. So, *^Zjv = ^7V) an d at the level of 
the integrable kernel: 

^ / x ^N{x l ]4> N {x 2 )-^ N {x l )(t) N {x2) 

K, N {x 1 ,x 2 ) = 

X\ — X2 

(t)_ N (x 1 )ip- N (x 2 ) - (j)- N (x 1 )^_ N (x 2 ) 
X\ - x 2 

= JC- N (x 2 ,xi) 

But, we can see on definition (Eqn. 1-7) that the correlators take a (—1)" sign if we 
revert the orientation of all the cycles. Thus: 

W„Oi, . . . , X n )- N = (-1)™ W n {Xi, X n ) N 

In particular, Wi must be odd in iV □ 
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